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ABSTRACT 

We use the BBGKY hierarchy equations to calculate, perturbatively, the lowest 
order nonlinear correction to the two point correlation and the pair velocity for 
Gaussian initial conditions in a critical density matter dominated cosmological model. 
We compare our results with the results obtained using the hydrodynamic equations 
which neglect pressure and we find that the two match, indicating that thare are no 
effects of multistreaming at this order of perturbation. We analytically study the 
effect of small scales on the large scales by calculating the nonlinear correction for a 
Dirac delta function initial two point correlation. We find that the induced two point 
correlation has a, behaviour at large separations. We have considered a class of 
initial conditions where the initial power spectrum at small k has the form A;" with 
< n < 3 and have numerically calculated the nonlinear correction to the two point 
correlation, its average over a sphere and the pair velocity over a large dynamical 
range. We find that at small separations the effect of the nonlinear term is to enhance 
the clustering whereas at intermediate scales it can act to either increase or decrease 
the clustering. At large scales we find a simple formula which gives a very good fit 
for the nonlinear correction in terms of the initial function. This formula explicitly 
exhibits the influence of small scales on large scales and because of this coupling the 
perturbative treatment breaks down at large scales much before one would expect it 
to if the nonlinearity were local in real space. We physically interpret this formula in 
terms of a simple diffusion process. We have also investigated the case n = and we 
find that it differs from the other cases in certain respects. We investigate a recently 
proposed scaling property of gravitational clustering and we find that the lowest order 
nonlinear terms cause deviations from the scaling relations which are strictly valid in 
the linear regime. The approximate validity of these relations in the nonlinear regime 
in N-body simulations cannot be understood at this order of evolution. 

Subject headings: Galaxies: Clustering - Large Scale Structure of the Universe, 
methods: analytical. 
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1. Introduction 

In the first paper of this series (Bharadwaj 1994), which we shall refer to as part I, we 
presented a formalism for perturbatively studying the growth of cosmological correlations using 
the BBGKY hierarchy. We also appplied this formalism to the particular case when the initial 
disturbance was Gaussian and could be described by the two point correlation function, all higher 
correlations being zero. We used the formalism to calculate, in a universe with Q = 1, the induced 
three point correlation function to the lowest order of perturbation theory. We found that the 
three point correlation function at any length scale depends on the two point correlation at all 
smaller scales i.e. the small scales influence the large scales. This disgrees with the commonly 
used 'hierarchical form' where one assumes that the three point correlation at any length scale 
depends only on the two point correlation at the same scales. 

The nonlinear evolution of the two point correlation function is driven by the three point 
correlation function through the tidal force of the third particle. In this paper we carry the 
perturbative scheme one step further and calculate the lowest order nonlinear correction to the 
two point correlation function. Before presenting this calculation we discuss some general issues 
about the scheme of calculation being used in this paper. We would also like to compare it with 
another scheme whose usage is more prevalent in the current literature. 

We consider a universe filled with coUisionless particles that are initially uniformly distributed 
and move according to Bubble's law. We study the statistical properties of disturbances in such 
an universe. To do this we construct an ensemble with each member corresponding to a different 
relization of the disturbance. We follow the evolution using Newtonian dynamics and use this to 
study the growth of correlation functions. The disturbances that are considered are such that 
initially they have a single streamed flow where all of the particles at any point have the same 
velocity and this velocity field is a smooth function of position. The initial mass density too is a 
smooth function . As the disturbance evolves the trajectories of particles intersect and the streams 
pass through one another interacting only through gravity. When this occurs the velocity and 
density are no longer smooth functions and they become singular and multi valued at points. This 
process is best illustrated by the Zeldovich approximation (Zeldovich 1970). Even though we start 
the evolution with a single streamed flow everywhere, multistreamed flows with particles with 
different velocites at the same point are inevitable. This fact forms the motive for the comparison 
of the two schemes of calculation given below. 

The scheme of calculation used in this paper is based on using the equation of motion 
of the individual particles to obtain equations for the evolution of ensemble averaged phase 
space distribution functions (i.e. the BBGKY hierarchy). We consider velocity moments of the 
distribution functions and evolve them perturbatively. This gives the evolution of correlation 
functions. These equations are valid in both the single stream and multi stream regimes. 

The scheme with which we would like to make the comparison evolves individual members 
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of the ensemble using the hydrodynamic equations i.e. the continuity equation, Euler's equation 
and Poisson's equation. Pressure and other higher moments of velocity, which are initially zero, 
are assumed to be zero at later times and this scheme may be described as the single stream 
approximation. These equations become invalid once multistreaming occurs. 

Although the two schemes are equivalent in the initial stages of the evolution when it is a 
single streamed flow, the single stream approximation breaks down as the evolution proceeds. 
Once this happens we expect the two schemes to give different results. The lowest order nonlinear 
correction to the two point correlation function has beenn calculated by Juszkewicz (1981), 
Vishniac ( 1983), and Makino, Sasaki and Suto (1992). Fry (1994) has calculated the higher order 
nonlinear corrections to the two point correlation function. All these calculations have been done 
using the single stream approximation. 

In section 2 we use the BBGKY hierarchy to calculate the lowest order nonlinear correction 
to the two point correlation function and compare it with the results obtained in the single stream 
approximation. In section 3 we calculate the lowest order nonlinear correction to the pair velocity. 

Unlike the linear evolution, the nonlinear evolution exhibits coupling of the various scales. 
This has been examined analytically by Suto and Sasaki (1991) and Makino et al. (1992). This 
issue has also been studied by Juszkewicz, Sonoda and Barrow (1984) and Hansel et. al. (1985), 
and for the particular case of the CDM initial condition it has been studied by Coles (1990), Jain 
and Bertschinger (1994), and Baugh and Efstathiou (1994). Most of the analysis has been in 
Fourier space. In this paper we address this question in real space. The coupling of the various 
scales in the nonlinear epoch has also been studied using N-body simulations by Mellot and 
Shandarin (1990), Little, Weinberg and Park (1991), and Evrard and Crone (1992). 

In section 4 we consider a situation where the initial two point correlation is zero beyond a 
certain separation (i.e. ${x) = for x > xq). In the linear evolution the two point correlation 
will remain zero for x > xq, but in the nonlinear evolution the two point correlation will develop 
because of the streaming of particles and the nonlocal nature of Newtonian gravity. This exhibits 
the influence of the small scales on the large scales. In section 5 we numerically investigate various 
cases where the initial two point function has a power law behavoiur x~'^ for large x. We look at 
three different quantities - the two point correlation, the average of the two point correlation over 
a sphere of radius x and the pair velocity. We look at both the spatial and temporal behaviour of 
these quantities. In part I we found that at large separations the induced three point correlation 
function has the 'hierarchical form' for a class of initial conditions, but there are initial conditions 
where, because of the effect of the small scales on the large scales, the induced three point 
correlation does not have the hierarchical form. We investigate if the nonlinear correction to 
the two point correlation function shows any qualitative differences for these two kinds of initial 
conditions. 

Hamilton et. al (1991) have suggested that in an f2 = 1 universe the evolution of the two point 
correlation funtion can be described by a simple scaling relation. More recently Nityananda and 
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Padmanabhan (1994) have examined the possible origin of this scaling relation. These arguments 
are based oninterpolation between the evolution in the linear regime and in the stable clustering 
regime. In section 5 we also check the validity of the proposed scaling relation in the weakly 
nonlinear regime. 

In section 6 we discuss a possible interpretation for an interesting result we obtain in the 
numerical investigations of section 5. This interpretation is in terms of a simple diffusion process. 
In section 7 we present the conclusions. 



2. Calculating the two point correlation. 



2.1. Notation and the equation governing the two point correlation. 



Throughout this paper we shall be using a comoving coordinate system. Instead of the cosmic 
time t, we shall use a parameter A defined by 

where S{t) is the cosmological scale factor. This parameter, which was also used in part I, was 
introduced by Doroshkevich et. al. (1980) and it make the calculations simpler. For an f2 = 1 
universe we have (part I) 

where p is the mean comoving density. We shall use Q = 1 throughout this paper. The momentum 
of a particle is defined by 

d 

We consider the initial density perturbation to be a Gaussian random field. This is completely 
described by the two point correlation function ${x). 

We also define a potential (j) such that in the linear regime 

a^\x') = l^Xct^(\x'-x'\)r.e' (4) 

where Aq is the initial value of the parameter A and e ~ (^6p/ p) is a small number in terms of which 
we are doing the perturbative expansion. 

We define f[p,\),c[x^,p^,x^,p^,\) and c?(l,2,3,A) as the reduced ensemble averaged one, 
two and three point distribution functions in the appropriate phase spaces. Notice that we have 
used (1,2,3) to indicate the point [x^,p^,x^,p^,x^,p^) in phase space. Sometimes we shall use 
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(1,2,3) to denote points in space. The difference will be clear from the context. The superscript 
indices a, b, c, etc. are used to indicate different positions. 

We shall be using the various velocity moments of these distribution functions. Some of them 
are defined below and the notation for the rest follow the same pattern. We have 

n = J f{l)dY {nm = p) , (5) 



Jplf{l)dV = 0, (6) 
where n is the number density of particles, 

n<{piy>i = J {piyf{l)dY-e', (7) 

n^^(x\x^) = J c{l,2)dYdY - e\ (8) 

<pl>2{x'x') = Jplc{l,2)dYdY^e' (9) 

<plpt>2{x\x') = Jplptc{l,2)dVdV-e\ (10) 

where a and b take values 1 and 2, and ^ is the two point correlation function. The subscript 
outside the angular bracket indicates which distribution function the moment refers to e.g. 
< >2 (1,2, A) is the first moment of the two point distribution function. The subscripts 
(/X, I/, a, etc) refer to Cartesian components and are to be summed over whenever they appear 
twice. Some more definitions are 

n^C (x\x^,x^') = J d{l,2,3)d^p^d^p^dY (11) 

n^<pl>3{x\x',x^) = Jpld{l,2,3)dVdVdV-e\ (12) 

where a takes values 1,2 and 3, and ^(1,2,3) is the three point correlation function. 

Also, we shall also use the following notation for the cartesian components of the inverse 
square force 

u < - x^„ 

Xf = (13) 

We next present the equation governing the perturbative evolution of the two point correlation 
function. This equation, which was derived in part I, is 



3 



f2-h- g^fl , (14) 
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where, 

d 



f, {1,2, \)=SGp— j C{l,2,iA)Xl^d^x\ (15) 
/2(l,2,A) = 25Grz— ^ / < >3 (1, 2, 3, A)Xf (16) 

Ox J. Ox,, J 



„1 „fc 

^ r " 

Here the position indices take values 1 and 2 and are to be summed when they appear twice. 



For an f2 = 1 universe this becomes 

24 5 24 d 
" J^dX^ + A^^ = " " 5A^^ ' • 

We separately consider the various terms on the right hand side of equation (18). We first 
consider equation (15). This depends on three point correlation function ( which has been 
calculated to order in part I. 

We reproduce the expression below. We have 

= ^(T)'^oE[3^(vV(a,a;)AvV(a,a;)) 



a=l 



(19) 



In the equation for the three point correlation function the following conventions are used 

A. the position indices, e.g. a, take values 1,2, and 3 corresponding to the corners of the triangle 
for which the three point correlation function is being evaluated. Also, a position index which 
appears twice or more should be summed over the allowed values. 

B. for a fixed value of the position index (e.g. a = 1), a'^ and a'2 are to be summed over the 
other two values (i.e. a'^ = 2,03 = 3 and a'^ = 3,03 = 2). This is to be done whenever such a 
combination of three position indices appear. 

In some of the equations for the other moments of the three point distribution function, if 
indicated, the summation convention A may not hold, but the convention B always holds. 

To calculate /2 and /s we have to first calculate the following quantities: < >3 (1,2,3, A) 
and d^d^d^ < p^^p^p^ >2 (1)2, A).. These calculations are are discussed in the following two 
subsections and the final result is presented in subsection 2.4. The overall strategy is the same as 
in part I. We take velocity moments of the BBGKY hierarchy and retain terms only up to order 
e^. As a result of this we obtain a set of ordinary differential equations in the parameter A. These 
equations have complicated spatial dependences. The spatial calculation is simplified by taking 
spatial derivatives (curl and divergence) and we then obtain a set of equations that can be easily 
solved simultaneously. 
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2.2. The triplet momentum. 



The triplet momentum < >3 (1,2,3, A) is defined as the first moment of the three-point 
distribution function c? (1, 2, 3, A). 

J p''^d{l,2,3,\)dYdYdY = < P],>3 (x\x\x^,\) . (20) 

It is a function of three positions 1,2 and 3, and the index 'a' in < >3, which indicates at which 
vertex of the triangle we are considering the momentum, can refer to any one of them. 

We want to calculate this quantity to order e^. This is the lowest order for which it has 
non-zero value for Gaussian initial conditions. 

The evolution of the triplet momentum is governed by the first moment of the third equation 
of the BBGKY hierarchy 



- SmGp j C (a",4)x^"(iV - SmGp j x (1, 2, 3, 4)X^"(i^ 

- SmGpi (a, «;) j i (4,4) X^^d^"" = . (21) 

To evaluate < >3 we separately consider both its curl and divergence with respect to cc" and 
use these to construct it. All the equations given below for the curl and divergence are valid only 
to order and in all of them the summation convention A does not hold. In all these equations 
the indices a, h and c refer to the three different corners of the triangle that we are considering i.e. 
a ^ h ^ c. 

The curl of equation (21) is 

of a 

+ ^ (y)'°^/^M^ dtV'cl> (a, 4)] = (22) 

where, 

F^"(l,2,3,A) = e^^,5;<p^>3 (23) 

and 

Gf{l, 2, 3, A) = e^^^d'^.dl < pIpI >3 . (24) 
and we have used the fact that for Gaussian initial conditions to order 



< pIpI >3 = < p1 >2 [a, aA <p1 >2 (a, aA (25) 
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to evaluate G^". 

The divergence of equation (21) is 



m 



Ao 



10 



I] ^5 (?) + ^5 (a) 



9=1 



where 



and, 



r{i,2,i,\) = di<pi >3, 

K''\l,2,^,\) = dldl<p1pl >3 



(26) 

(27) 
(28) 



+25^5^ W (a, «;) 5^5^ W (a, 4) • 



(29) 



We have used equation (26) to evaluate if" and equation (19) for the three point correlation 
function. 



Next we consider the second moment of the three point distribution function 

^ < p>J; >3 J- 



^5,. <!'Xf5>3(1.2.3.A) 



Sm'^G 



+ Cvl) f (c) d [c , 3] Xxt^d'x''d''p 



Sw?G 



+ Cvl) c ( c, ci j c ( C2, 4j X^Wd^'p 



-4c 73„4 7l2„ 



/ (O' + •^f.P;:) e (1, 2, 3, 4)X,4'^ci^a;4cii2^ = . 



(30) 



We use this to get an equation for 



9 riab 



abc , 6m 



Sm^^e^^. [5W</.(a,6)5;vV(a,c) 
5;vV(a,&)5."VV(a,c)] =0. 



where 



We have used the fact that for Gaussian initial conditions to order 

<pX?''>3 = <P>l>2{a,b)<pl>2{a,c) 



(31) 
(32) 



<P>i >2 {a,h)<pl >2 (a,c) 



(33) 
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to evaluate H'^'^ when d = a oi d = b. 

By taking divergences of equation (30) we obtain 



^ l^a^c + ^ (j'' + J^) + m^^ [5 (a) + 5 (6)] = (34) 



where 

rabc Qaabac ^ ^ci 



d^dtd'. < pXpI >3 • (35) 



Finally we have the third moment of the equation for the three point distribution function 

^ < P>tpl >3 < P>Up', >3 (1, 2, 3, A) (36) 

[S;>tpl + ^tPlPl + f (e) d (e" , 4) X^^d'x'd''p 

[^ZpIpI + Cp>^^ + ^ (e> ^i) ^ (4,4) X^^d'x'd'^p 



Sm?G 



Sw?G 



Sw?G 



j {^ZpIpI + Cp>^^ + e (1> 2, 3, 4)X^^d'x'd''p = . 



This can be used to obtain the equation for H^^'^ 

^ Tjabc I ( r^ab , r^ac 

+ lOm^ ep^^ [5^VV (a, &) ^^VV («, c) 

+ 5;vV(a,&)5."VV(a,c)] =0 (37) 
To obtain this equation we have used the fact that for Gaussian initial conditions to order 

< P^Plplp'-y >3 = < P^pI >2< pIp"^ >3 +permutations . (38) 
Taking divergence of equation (37) we obtain 

^afcc ^ O /^afc ^ ^fcc ^ ^caA + 3 ^ mM -f^ U (p) = . (39) 



d\ A2 V y ' ^ I ^12 

The equations (22), (31) and (37) can be simultaneously solved to obtain 

m /A^°^ 



Fl = f (^^j6/3^.[5;vV(a,a;)5,"VV(a,4)] (40) 
Of = m2(^^^6^^.[5;vV(a,a;)5."VV(a,a;)] (41) 
Hf^ = 2m' (^] ep^, [S^VV («, a,) d^V'cj, (a, 4)] . (42) 
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Simultaneously solving equations (26), (34) and (39) we have 



14 I A9 



[2g{a) + g{b) + g{c) 



K 



ah 



m 



r abc 



■^\:^^\[2g[a) + 2g[h) + g[c) 



Using these we obtain the triplet momentum as 



\10 r-i 

< pl >3 (1, 2, 3, A) = [-5;VV [a, a^) [a, 4 



-d;[d:V'(l^[a,a,)d:V'(l^[a,a, 

^5;vv vv + ^5;5^vv 5^vv 

^d^dt'cl) (a, a;) 5^VV (^1,4) + ^d^dt'd^'^cl) (a, 4) S^S^^V (4, 4 
^ / ^r^.' (^.'VV (4, 4) W (4, 4)) ^'-^ 



567r 

and we also obtain for a ^ b 

™2 \10 

dt < pIpI >3 (1, 2, 3, A) = — ^ [75;VV (a, 4) VV (a, 4 
+ 25;(5^VV(a,4)5^VV(a,4)) +205;vV(a,6)VV(&,c) 
+ 105;VV (a, c) VV (&, c) + 145;5^VV (a, 6) 5^VV (&, c) 
+ 145;5^</. (a, 6) V (&, c) + 75^5^ W (a, c) 5^ V ^ (6, c) 
+ ldldlci> (a, c) 5^VV (&, c) + 85;5^5^</. (a, h) dld^c^ (6, c) 



(43) 

(44) 
(45) 



(46) 



(47) 



2.3. The third moment of the two point distribution function. 



In this section the position indices takes the values 1 and 2, and they are to be summed 
whenever they appear twice or more. 

The third moment of the two point distribution function is governed by the equation 
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Sw?G 



Sw?G 



J {^ZpIpI + ^tPlPl + d (1, 2, ^)X^^^d^xH^p = . 



If we take divergence with respect to all the three sets of free indices we have 

^d;dtd^^ < pIpIpI >2 +^d;dtd^^d^ < p^pIpIp'!, >2 

+ U-KSmGp [< pIpI >i dldti{l, 2, A) + d^dt < p^pl >2 (1, 2, A)] 

- ^SmCpd^dldl j X^^ < pIpI >3 (1, 2, 3, \)d^x^ = . (48) 

For Gaussian initial conditions the terms in this equation have non-zero values only at order 
and higher. To order there are two unknown functions in the equation i.e. d^d^d^ < p'^p^p% >2 
and < p^p" >2- We need one more equation to self-consistently solve this equations. This equation 
is obtained from the second moment of the two point distribution function 



^ < P»2(1,2,A)+-— <p»:>2(l,2,A) 

- ^ / {C,pI + Cpi) f{c)c{c, 3) X'^^d\'d'p 

- '^J{6llpt + 6tpi)d{l,2,3)X'^^d'x'd'p = 0. (49) 
From equation (49) we get 

^d;d^ < pIpI >2 +^d;dtd^^ < p^pIpI >2 

+ ^Wdt < pIpIpI >2 -2SmGpdldl \ <pl>^ Xfd'x' = . (50) 
Simultaneously solving equations (48) and (50) we obtain 

d^^dldi < pIpIpI >2= ^d;,d:d: < pIpIpI >2 +^d;,dtdtd^ < pIpIpIp", >2 

+ 26A < >^ Wai, 2, A) - ^;^dld: j <pl>, X^'^d'x' 

- ^d^dtd^^ J < pIpI >3 Xl^d^x^ (51) 

All the terms on the left hand side of equation (51) are known to order e^. Writing it in terms of 
^ we have 



d^dlK < pIpIpI >2= ( ^ 



82 



^ VV(1,2)VV(0) 



- 12 - 



240 .^.^ ,xv74^n o^ , 816 ^i„4 



— V> (1, 2) V> (1, 2) + — 5X</. (1, 2) SX*/- (1> 2) 

965;.5>V(i,2)5;.5^vV(i,2) + — 5;.5XvV(i,2)5;.5^5i</.(i,2) 
485;.5>V(i,2)5;.5^</.(i,2) + y5;.5X5^</'(i,2)5;.5^5i5^</.(i,^ 

/ Xf (365>V (1, 3) VV (2, 3) + f S^S^VV (1, 3) 5^VV (2, 3) 



45 



-5i5^</. (1, 3) 5^VV (2, 3) + 9dtdid$cl> (1, 3) 5^5^ V^</. (2, 3) d'x 



(52) 



2.4. The two point correlation. 



Using the results derived in the two previous sections we can calculate /2 and /s. We first 
consider only the A dependence of ^. Equation (18) may be written as 



5A3 

This has a solution 



24 5 24 
A2 5A^+ 



_£ 
A 



/2(Ao)-/3(Ao) + ^/i(Ao) 

^0 



ai,2,A) 



504 A8 



/2(Ao)-/3(Ao) + ^/i(Ao) 

^0 



(53) 



(54) 



Using equations (15), (16) and (17) we calculate the spatial dependence of the right hand side. 
This gives us 

1 /a1°\ r 

^(2) (1, 2, A) = — f ^ j [-32VV (1, 2) VV (0) - 14VV (1, 2) VV (0) 
+ ^VV (1, 2) VV (1, 2) + 3605X'/' (1> 2) dlV^ (1, 2) 

+ ^5;.5^vV(i,2)5;.5>V(i,2) + i445;.5XvV(i,2)5;.5^5i</.(i,2) 

1 44 

+ i265;.5^vV(i,2)5;.5^</.(i,2) + — 5;.5^5i5^</.(i,2)5;.5X5^</'(i,2) 

- ^/x^^g5ivV(l,3)5^5^VV(2,3) + ^5i5>V(l,3)5X5^VV(2,3) 



15 



2 5i5>(l, 3) 5^5^5^VV (2, 3) + Zdldldlcj) (1, 3) V^</. (2, 3) 

21 



- 15V> (1, 3) 5^V> (2, 3) - y S^VV (1, 3) didlV'cl> (2, 3) ) d'x 
This is the two point correlation function at order e^. 



(55) 



We do a Fourier transform of (55) and compare it with the result obtained in the single 
stream approximation (Makino et. al 1992) and find that the two match. 
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The algebra involved in deriving equation (55) was checked using the mathematical package 
MATHEMATICA and the Fourier transform was done using this package. 

The calculation presented here, which is based on the equations of the BBGKY hierarchy, 
has the effects of multistreaming, if any, at the lowest order of non-linearity i.e. e^. This matches 
with the results obtained in the single stream approximation which does not take into account any 
effect of multistreaming. Hence we find that there are no effects of multistreaming at this order of 
non-linearity. 

Given the initial two point correlation function ^^^\x, A), to evaluate the non-linear correction 
at order we solve the equation 

V2(VV(a;)) = f^i(a;,Ao) (56) 

to obtain V^^(a;) and then solve for 4>{x) and use these in equation (55). These calculations 
are simplified a lot if we use the fact that ^(cc, Aq) is spherically symmetric. Equation (56) does 
not uniquely determine the functions V^^(a;) and 4>{x). If V^^(a;) is a solution of equation 
(56) then so is V^^(a;) + C where C is some constant. Under this transformation we also have 
^(cc) ^(2;) + ^f-- If we consider all the terms in equation (55) that are affected by this 

- 32VV (1, 2) VV (0) + 126dldlvU (1, 2) dXcj, (1, 2) 

- ^Jx'^'dldicl>{l,3)dXdlvU{2,3)d'x' (57) 

we find that the correction to the two point correlation function is unchanged by such 
transformations and independent of C. Similarly, we can add a constant to 4>{x) but this obviously 
does not have any affect as only derivatives of 4>{x) appear in equation (55). Hereafter we shall 
use the boundary condition that V^^(a;) should vanish as x goes to infinity to fix the constant C. 

At this stage we should point out that it such a choice of C is not always convenient. For 
example if the initial two point correlation is such that the power spectrum has the form P{k) oc A;" 
with n < —1 at small k, then the boundary condition stated above implies that 

poo 

VV(0) = -47r / P{k)dk (58) 
Jo 

is infinite. Although V^^(O) is infinite and this quantity appears in the non-linear correction to 
the two point correlation function ^^^\x,t), we may still get a finite ^^^\x,t) for certain initial 
conditions. This is because now is the difference of two infinite large quantities which cancel 
out to give a finite result. The same problem is encountered if one does the analysis in Fourier 
space where for —l<n>—3 the correction to the power spectrum is a finite quantity which 
is the difference of two divergent integrals (Vishniac 1983). In real space this situation is easily 
handled by changing the boundary condition used to calculate V^^(a;). If we use the boundary 
condition V^^(O) = to fix the constant C than the situation discussed above does not occur and 
it is possible to calculate ^^^\x,t) solely in terms of finite quantities. Here we shall only deal with 
situations where the former boundary condition ( limj;_>oo V^^(a;) = 0) can be applied. 
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3. The pair velocity 

We next calculate the first moment of the two point distribution function, < >2 (1,2), 
to order . This is a function of two positions and the index a indicates at which of the two 
positions the momentum is being considered. We use the pair continuity equation 



^ i[l,2,\)+^-^<pl >2 (1,2,A) = 0. 



m 



to obtain 



where 



(59) 



(60) 



a a 

2/^ = 



end d refers to the complement of a (e.g. if a=l, d =2). Using the two point correlation at order 
calculated in the previous section we get at order 



<PI >2 
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A quantity related to the first moment of the two point distribution function is the pair current 

<pI>2- <pI >2 



(61) 



m 
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whose divergence gives the rate at which the correlation at any separation is growing. We use 
this to calculate the pair velocity which is the ensemble average of the relative peculiar velocity 
between any two particles at a comoving separation x at time t (or A). 

In terms of the pair current this is 

.,(.,0 = 5^ 
4. Effect of small scales on large scales. 



wiix^X), 



(63) 



In this section we study how the small scales affect the large scales. We consider a situation 
where initially ^{x) = for x > xq. In the linear regime ^(cc) will continue to have value zero for 
X > xq. This will not be true in the nonlinear regime and there will be nonzero correlation for 
X > Xq due to the streaming of particles and the non-local nature of gravity. In this section we 
want to find out, to the lowesst order of nonlinearity, the nature of the induced correlations. 

We first consider an extreme case of the above situation where the initial two point correlation 
is a Dirac delta function and ^^{x) = for x > 0. We initially have 

^^'\x, Ao) = ^AoVV(2;) = \\oA6\x). (64) 
Solving for the potential we get 

cl>{x) = -^x (65) 
and using this in equation (55) we obtain for x > 0. 

This shows the influence of the small scales on the large scales. 

We next write equation (66) in a form that can be compared with a general case where the 
initial two point correlation function has compact support. This is done by introducing a quantity 
M(A) which is related to the integral of ^(■^^(a;, A) over the whole region x < xq where it is non-zero 
i.e. 

M(A) = J ^^^\x,\)x^dx. (67) 
We can then write equation (66) for x ^ xq as 

^(2) (^, A) = 3.264^^. (68) 
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The general case where the initial two point correlation has compact support cannot be 
treated analytically because the integrals in equation (55) cannot be done analytically. 

To check whether the conclusions drawn from the case where the initial condition is a delta 
function can be generalized we have considered a situation where 

V^(l){x) = 16 - 12a; + for x < 2 (69) 

end 

VV(2;) = for X > 2 . (70) 

This corresponds to the self- convolution of a spheres of unit radius. The corresponding power 
spectrum has the form 



P{k) (X 



sin(A;) — k cos(A;) 



2 



(71) 



A;3 

Doing the integrals in equation ( 55) numerically we find that for large x i.e. in the interval 
30 < cc < 40, the induced two point correlation can be fitted by the form 

^(2) (^, A) = 3.277^^. (72) 

X 

For larger values of x the integrals could not be reliably evaluated numerically as the numbers 
involved are extremely small. The initial and induced two point correlation are shown in figure 1. 

We propose that the form (68) holds for the two point correlation induced ar large distances in 
all cases where the initial two point correlation has compact support and a nonzero total integral. 

We next consider the pair current. The initial pair current for the delta function initial 
condition is 

3,{x,\) = l^M(\) (73) 

This current is divergence free everywhere except at the origin. Since it does not change the 
correlation at any nonzero separation it may be thought of as matter flowing into the overdense 
regions from infinity and flowing out of the overdense regions to infinity (Zeldovich and Novikov 
1983, section 10.4). 

If we try to calculate the nonlinear correction to the pair current for the delta function initial 
condition we get a divergent answer. This is because the integral /q°° i^^\y)y^dy diverges . 

For the other initial condition considered above (equation 71) we find numerically that the 
induced pair current (figure 2) has a x~^ behaviour at large x. This is the leading part of the 
induced pair current and it too represents the flow of matter from infinity. In addition to this 
divergence-free part, the pair current also has a part that has a x~^ behaviour at large x and it 
corresponds to a local redistribution of matter. It is the latter that gives rise to the x~^ correlation 
at large x. 

For large separations (cc > 2) the pair velocity has the same spatial dependence as the pair 
current. 
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Peebles (1980) and Zeldovich and Novikov (1983) (and references therein) have discussed the 
effect of small scales on large scales. They have considered an initial power spectrum which has 
the value zero at small k (large length scales). They find that the rearrangement of matter on 
the small scales gives rise to a power spectrum of the form P{k) a A;^ at small k or large scales. 
N-body simulations by Mellot and Shandarin (1990) find that there is indeed a A;^ power spectrum 
generated in the initial stages of the evolution but this changes over to a k^ form by the time 
the small scales saturate. The initial condition considered in these calculations corresponds to a 
situation where the universe can be divided into small cells and the matter content in all the cells 
is exactly equal. Also the matter in each cell is distributed such that its center of mass coincides 
with the center of the cell and the total momentun of each cell is zero. They have investigated 
the effect of redistributing the matter inside each cell, conserving the mass and momentum of the 
individual cells. Matter and momentum conservation are of course also built into the scheme used 
in this paper, but the analysis presented here differs in that the initial conditions used here are 
such that the cells all have different masses and the mass fluctuations of cells that are at large 
separations are uncorrelated. These are not the same initial conditions that give rise to the A;^ 
spectrum. We have studied the nature of the correlations that develops between the fluctuation in 
the mass of cells that are at a large separation. Note that this is a quantity in real space whereas 
the A;^ spectrum being analytic in the limit A; ^ does not imply any long range power law 
correlation in real space (e.g. P{k) = k^e~^^ has a Fourier transorm V^e~^^/^. Hence the results 
given here are not directly related to those of Peebles (1980). 

5. Numerical Investigations. 

The expressions for the contribution to the two point correlation and the pair velocity at order 
are rather complicated. To understand them better we proceed to evaluate them numerically 
for different initial conditions. 

5.1. The spatial behaviour of the contributions. 

First we consider initial two point correlations such that the corresponding power spectrum 
has the form P{k) = k'^e~^. We consider cases where n takes the values n = .5,1,1.5,2 and 
3. These initial conditions have just one length scale which is introduced by the exponential 
cut-off which becomes effective for A; > 1. In all these cases the correlation function has a 
power law behaviour x~'^ for large separations and the corresponding values of the index 7 are 
7 = 3.5,4,4.5,6 and 6. We also consider the case P{k) = ke~^^ which has a Gaussian cut-off for 
the power spectrum at large k instead of the exponential cut-off used in all the other cases. 

Figures 3 shows the function i^^\x) for two of the cases. This function is defined so that at 
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any instant the two point correlationn at order is 



Figure 3 also shows the quantity ^(^)(a;). This when multiplied by ^^^j gives the ~ 
contribution to ${x,t) which is defined as 

i{x,t)=\ r i{y,t)yHy. (75) 
Jo 

This is the average of the two point correlation function over a sphere of radius x . All of the 
cases considered here satisfy the relation $^^\x) oc $^^\x) for large x. This is crucial in deciding 
the behaviour of the induced three point correlation function ^(1,2,3) at large separations and for 
this class of initial conditions we have ^(1,2,3) ~ ^(1,2)^(2,3) + permutations. 

Figures 4 shows the function v^^\x) which is related to the radial component of the pair 
velocity at order at any instant by the relation 

For all these cases the initial pair velocity too has a power law behaviour x~^ at large 
separations and we have /3 = 7 — 1 as expected. 

It should be pointed out that the initial conditions when n = 3 differs from the other cases. 
For all the other cases $^^\x) crosses zero only once, it is positive for small x and goes over to a 
negative value at lerge x where it has the power law form. The quantities ^^^\x) and v'' ^\x) do 
not change sign , the former is positive and the latter is negative. In the case when n = 3, $^^\x) 
crosses zero twice, and $^^\x) and v^^\x) cross zero once. Thus in this case (tc = 3) at large x the 
signs of all the quantities are opposite to the signs in the other cases. 

We have calculated the contribution to all of the above mentioned quantities in the renge 
< cc < 40. 



We first discuss the large separation behaviour of $^^\x) which is defined such that the 
contribution to the two point correlation is ^ ^ j This function is shown in figure 5. 

For all the cases we find that at large separations $^^\x) has a power law form $^^\x) ~ x~'^ with 
7/ = 7 — 2 i.e. ^(2)(a;) ~ V2^(^)(a;). Motivated by this we investigated whether there is any simple 
relation between ^^^\x) and ^^^\x) which holds for all the cases. We looked at the ratio 



V2</)(0)V6</)(a;) 

at large x for the different cases and we find that the value of Z is nearly the same 
( — .048 < Z < —.049) for all the cases. This relation can also be expressed as 



(77) 



r 7"°° 1 

^^'\x) = R / ^^'\y)ydy V'^^'\x) (78) 
Jo 
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with R = —4Z ~ .194. In terms of the power spectrum this may be written as 

R 

2^ 



P2{k) 



f 

Jo 



Pi[k )dk 



(79) 



where Pi{k) is the initial power spectrum and P2{k) is the correction to the power spectrum at e^. 
The fact that /q°° P{k)dk > tells us that P2{k) < 0. In real space we can say that (P'\x) has the 
same sign as (y\x). 

Based on this numerical evidence we make a hypothesis that equation (78) holds for all initial 
conditions where the initial power spectrum index n satisfies < tc. We have not considered any 
cases where n is negative. For such cases in < —1) the integral in equation (79) does not converge 
and we do not expect this relation to hold. 

One of the factors we believe to be responsible for the spread in the values of R (or Z) is that 
i^^\x) goes to the power law form asymptotically and for all the cases we have not been able 
to calculate i^^\x) to equally large separations. As of now we have no rigorous derivation for 
equation (78), but we give a heuristic interpretation in terms of a diffusion process in section 6. 

Makino et. al. (1992) have analytically calculated P2{k) for various power law initial 
conditions, of which the case where 

. k 



mi 



kr 



for k < kr 



(80) 



and zero elsewhere is of interest to us. In the limit k— > their analytic expression for P2{k) 
reduces to 



This can be written as 



61 / 1 
315 V2^ 
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Pi[k )dk 



k^Pi{k) 



which on comparison with equation (79) gives us i2 = ^ ot Z = — ^ i'ih) ~ ~-0484. 



(81) 
(82) 



We see that this matches equation (79) and this serves as a test of our hypothesis. 

At small X the behaviour of $^^\x), as shown in figure 5 is rather complicated and is difficult 
to generalize. We find that it starts with a positive value at cc = 0, which falls fast, becomes 
negative and oscillates around zero a few times before going over to the power law fowm. 

We next consider the behaviour of ^(^x,t). At large separations this has a behaviour quite 
similar to ^^^\x) as shown in figures 5 and we find that 



r poo 1 1 

e'){x) = R / ^^'\y)ydy -V\x^W{x)) . 
.Jo J X 



(83) 



We see that at large separations ^(^)(a;) has the same sign as ^(^)(a;). The small x behaviour of 
^(2)(a;) is similar to ^^^\x). 
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The behaviour of v^'^\x), which is defined such that the contribution to the pair velocity is 
5(f)^5(f)S(2)(a;), is shown in figure 6. For large values of x the behaviour can be described by 

[poo 1 

v^^\x) = R / $^^\y)ydy V^v^^\x) . (84) 

We see that v^'^\x) too has the same sign as v^^\x) at large x. At small x the function v^'^\x) 
starts from zero and rises fast and then falls off and changes sign a few times before going over to 
the power law form. 

We see that at very small scales the contribution acts to increase the correlation, whereas 
at intermediate scales it can act to increase oe decrease the correlations.. 

We next consider a case where the initial power spectrum has the form Pi{k) = e~^. The 
two point correlation function has the form $^^\x) ~ x~^ for large x but $^^\x) does not have 
the same behaviour and we have ^^^\x) ~ x~^ instead (figure 7). Because of this, at large 
separations the three point correlation function exhibits a behaviour which is quite different from 
the one exhibited by the cases considered previously and in this case the three point correlation, 
as discussed in part I, does not have the 'hierarchical form'. We find that the large x behaviour 
of ^^^\x) too is somewhat different as compared to the previous cases. Although we find that 
^^^\x) oa V^^(-^)(a;) the factor relating the two is diferent from that found for the previous cases 
and if we fit a formula like equation (79) we get R ~ .496 instead of .194. 

We next consider the large x behaviour of ^(^)(a;). We see that it has a behaviour of the form 
x~^ as compared with the initial function ^^^^{x) which has the form x-^ i.e. a difference of three 
in the power law index. This is different from the previous cases where there was a difference of 
two between the index for the linear function and the nonlinear correction as seen in equation 
(83). This can be easily understood by noting that if we try to relate ^(^)(a;) with ^^^^{x) using an 
expression like equation (83) we find that the fact that ^^^^{x) ~ x~^ implies that the right hand 
side is zero. We then deduce that, in addition to a part that behaves as x~^, ^^^^{x) has part that 
behaves as cc"^. At large x the value of ^^^^{x) is determined solely by the former term as the 
latter fall off much faster, but the behaviour of the nonlinear correction ^^^\x) is determined by 
the latter as the first term does not contribute. We find that equation (83) gives a good fit for 
R ~ .496 which is consistent with the fit for ^^^\x). 

At large x the behaviour of v^'^\x) is similar to the behaviour of ^(^)(a;). 

At small X both $^^\x) and $^'^\x) start of with positive values whichh fall off fast. At 
intermediate values $^^\x) changes sign twice and goes over to the the power law form at large x 
whereas ^(^)(a;) changes sign only once and hence the two quantities have opposite signs at large 
X. The behaviour of v^'^\x) is similar to the behaviour of $^^\x) axcept that it starts from the 
value zero. 



5.2. The temporal behaviour 
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Here we would like to investigate the evolution of the correlation function and the pair 
velocity. It is generally believed that the linear results should hold on some legthscale until 
the density contrast averaged over that length scale is of the order of unity. Having calculated 
the lowest order nonlinear term we can see when this becomes of the order of the linear term. 
This would be a different criterion for determining when the linear results would no longer be 
applicable. Here we wish to compare these two criteria and investigate whether they are the same. 

Before proceeding further we should remind the reader that we are working in the continuum 
(or fluid ) limit and the initial conditions are such that the perturbative treatment is valid at 
all length scales. As a result of this there is a growth of clustering even at the smallest scales. 
Although in reality there are bound particles at some extremely small scales these are not taken 
into account. 

We only discuss the case with Pi{k) <x ke~^ here. We find that the other cases considered 
have a similar behaviour. 

The smaller scales go nonlinear first. We first consider the evolution of ${0,t) which is the 
mean square density fluctuation (figure 8). We find that the nonlinear correction enhances the 
growth of the correlation. We also find that the linear term is equal to the correction when ${0,t) 
is of the order of unity. This is as expected and here it happens at S{t) = .121. At S{t) = .053 
the correction is one tenth of the linear term and we may expect that the other higher order 
corrections not considered here will not contribute before this epoch. 

Next we consider the separation x = .1. Figure 9 shows the evolution of the various quantities 
of interest. There is no qualitative difference here with x = except that at cc = the pair 
velocity is zero. When x = 1 (figure 10), there is a qualitative difference. Previously the effect 
of the nonlinear terms was to increase the correlation and now it tends to decrease it, but here 
too ^(^)[x,t) and ^(^)(a;,f) are equal when they both are of order unity. At very large separations 
e.g. X = 20 (figure 11) the behaviour is quite different. We find that ^(^)[x,t) and ^(^)(a;,f) are 
equal when ${x,t) ~ .01. The correlation function and the pair velocity show a similar behaviour 
too. Hence it appears that at large scales the linear theory is breaking down much before one 
would expect it to. It should be noted that this happens at a very large value of S{t) (~ 10) 
and the perturbative approach has broken down much earlier on the small scales. The correction 
^(^)(a;,f) = R ^^^\y,t)ydy^ V^^^^\x,t) is nonlocal and it has a relatively large contribution 
from the small scales. In this expression the contribution from the small scale keeps on growing, 
whereas in reality the clustering saturates at the small scales because of virialization. Peebles 
(1980) has argued that once virialized objects have formed on some small scales, those small scales 
have no influence on the evolution of the large scales. Taking this into account we see that at large 
scales the lowest order nonlinear correction overestimates the contribution from the small scales in 
the epoch when the small scales have become highly nonlinear and have formed structures which 
are in virial equilibrium. We are still left with the fact that any particular small scale affects the 
evolution of the large scales until the small scale in question gets virialized. Just how significant 
this effect is has to be estimated by some independent means (i.e. not lowest order nonlinear 



- 22 - 



perturbation). One possible method of doing this is to compare the perturbative results with 
N-body simulations. Little, Weinberg and Park (1991), and Evrard and Crone (1992) have used 
N-body simulations to study how various scales affect one another in the non-linear epoch for an 
initial power spectrum of the form P{k) a . They find that the large scales have a strong 
influence on the evolution of the small scales but the small scales have very little influence on the 
evolution of the large scales. This is not in direct conflict with the findings of this paper because 
they are not claimed to be valid for the n = —1 type spectra. 

In a more recent paper Jain and Bertschinger (1994) have studied how the various scales 
influence one another in the nonlinear epoch for the CDM initial conditions where the linear power 
spectra has Pi{k) oc k for small k. They have calculated the lowest order nonlinear correction to 
the power spectrum P2{k) and they compare it to the results from N-body simulations. They 
follow the evolution until the density contrast on the scale of 8h~^Mpc is of the order of unity and 
they find that at small k the nonlinear correction to the power spectrum is negative and is much 
smaller than the linear power spectrum {i.e.P2{k,t) <C Pi{k,t), which is in good agreement with 
their N-body simulations. We have considered the CDM spectrum with the normalisation used by 
Jain and Bertschinger (1994) and used it in equation (79) to calculate the lowest order nonlinear 
correction to the power spectrum. The results, not given in detail here, are consisitent with those 
obtained by Jain and Bertschinger(1994) {i.e.P2{k,t) <C Pi{k,t) . For the standard CDM initial 
conditions with the normalization ag = 1, the condition ^^^\x,t) ~ ^(^)(a;,f) <C 1 at large x will 
occur only if we extend the evolution into the future. 

In section 6 we discuss a possible physical picture for the influence of the small scales on the 
large scales. We also discuss a possible means of independently estimating how significant this 
effect will be for general initial conditions. 

5.3. Scaling relations. 

Hamilton et. al. (1991) have suggested that in an f2 = 1 universe the evolution of the two 
point correlation funtion can be described by a simple universal relation whose exact form they 
have obtained by fitting N-body simulations. More recently Nityananda and Padmanabhan (1994) 
have examined the possible origin of this universal scaling relation. The scaling relation can be 
based on the conjecture that the dimensionless pair velocity h[a,x) = — depends on S{t) and 

X through ${x,t) alone. This conjecture is valid in the linear regime. Here we perturbatively test 
this conjecture at the lowest order of nonlinearity. 

We look at the behaviour of h(^x,t) and ${x,t) in some range where we are reasonably sure 
that the perturbative approach gives a good description of the clustering. Since the smallest scales 
go nonlinear first, the criterion is based on the properties at cc = and we restrict our analysis to 
the epoch when ^^^\o,t) < .l^^^\0,t). B ecause of this conservative criterion we find that h(^x,t) 
differs very little from the linear value of |^(a;,f). To look at the nature of this small deviation 
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we consider the ratio iifi!). Figure 12 shows this ratio as a function of ${x,t) for different S{t) 

for the initial condition Pi{k) oa ke~^. We find that for a fixed value of ^{x,t) the ratio has a 
spread of values. Although this spread is not large (5 percent), it is comparable to the largest 
deviation of mean square density from the linear prediction (ten percent). We also find that this 
ratio has a systematic behaviour which can be understood by looking at the same quantity as a 
function of x for different S{t) (figure 13). We see that there are points where the corrections to 
the pair velocity v^'^\x) and also to ^(^)(a;) are both zero and h(^x,t) and ^{x,t) there (figure 12) 
continues to follow the linear evolution. If we consider the first such point, then for smaller values 
of X (larger values of ^(^x,t)) the deviation is positive with respect to the linear value. For larger 
values of x (smaller values of ^(^x,t)) the deviation is negative with respect to the linear value. 
At large x (beyond the second zero crossing of the correction to the pair velocity) the deviation 
again becomes positive, but this is not seen in figure 12 because the whole range of x gets mapped 
to a very small range of ${x,t) between and some very small value. Based on this we draw the 
conclusion that we cannot express h(^x,t) as a function of ^(a;,f) alone. 

We have carried out this excercise for all the different initial conditions discusssed earlier and 
they too exhibit a simiiar behaviour. 

As was discussed earlier, equations (83) and (84) give a good description of ^(^)(a;) and v^'^\x) 
at large x for a large set of initial conditions of the form ^^^\x) oa x~'^ . Using these and the fact 
that ^^^^ = — i^-ll^ we have 

X 3 

h{x,t) = 2S{tf 

and 

i{x,t) = s{tf 

This also show that h(^x,t) is not a function of ${x,t) alone. 

6. An interpretation based on diffusion. 

In the previous section we saw that for a certain class of initial conditions equation (78) gives 
a good fit for 

in terms of ^^^\x). In this section we provide a possible interpretation for 
this equation. The cosmic energy equation (Peebles 1980 and reference therein), which is the 
second moment of the first equation of the BBGKY hierarchy, allows us to relate the integral that 
appears in equation (78) to the mean square momentum. In the notation used in this paper the 
cosmic energy equation can be written as 

d d f°° 

g^<{py >iW = 47rSGpm^—J i{x,X)xdx (87) 



2S{tf^{^-l)R 
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where x =\ — x'^ \. 

Instead of looking at the evolution of < {p^Y >, we consider the motion of the particles as a 
function of the growing mode which in this case is the scale factor S{t). We define a velocity 

db mdb 

and < > {S{t)) is the mean square of this velocity. In terms of this the cosmic energy equation 
is 

r„, o 1 3S d 

'dS 

which in the linear regime gives us 



r 1 n r°° 

[s'<u'>] =Y-dsJo 



1 P'^ 

< >= ^('\x)xdx, (90) 

i.e. < > does not change. Note that in this equation we have not explicitly shown the 
superscript (1) over u to indicate that it the linear part of u. Henceforth we shall use u for the 
linear part of the same quantity. Using equation (90) we can write equation (78) as 

^('\x,t) = RS\t) <u'> (1^) ' V'^('\x) (91) 

which looks like the solution of a diffusion equation. It is in light of this that we interpret equation 
(78). 

Consider a particular realization of the Gasussian density fluctuation field. We consider two 
points x^ and x^ where the density fluctuations are A[x^,t) and A[x'^,t) respectively. The A[x,t) 
of the fluid element at any point x grows according to linear theory and the fluid element moves 
according to some random velocity which is assumed uncorrelated to the velocity at any other 
point. We then have 

S^(t) 

+ S{t)Uf,df,A{x,to) + —^Uf,u^df,d,.A{x,to) . (92) 
Using this and taking an ensemble average we have 

ax'-x',t) = <A{x\t)A{x',t)>=(J^ya^'-x',to) 

+ ^S'{t)<u'> ^-^^\'a^'-x\to) (93) 

where 

< u^u^ >=^J^<v? > . (94) 
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The first term on the right hand side corresponds to the linear growth and the second term is the 
lowest nonlinear correction. We see that by this process we have a nonlinear correction which 
matches with equation (91) except for the numerical factor R. This difference can be attributed 
to the fact that we have kept the velocity li^ constant and this is only valid in the linear regime. 
In the nonlinear evolution this is not true. Based on equation (93) we suggest that the simple 
diffusion process described above gives a description of the evolution of the two point correlation at 
the lowest order of nonlinearity for large separations. Thus, when we are looking at the correlation 
at a large separation x, we can consider the initial perturbation field to be growing linearly, but 
there is one more effect to be taken into account. This effect is the local rearrangement of matter 
on small scales and this may be thought of as being random. The local rearrangement is on a 
lengthscale L ~ S{t)V < v? > and we expect this picture to be valid when L <^ x. When the 
mean displacement becomes comparable to the separation {i.e.L ~ x) we can longer treat the 
displacements as random and this picture is not valid. 

We see that X2(f), the mean square displacement of the particles from their initial positions 
is the quantity that determones the influence of the small scales on the large scales. It is also 
this quantity that gets overestimated in the later epochs in the perturbative result. To make an 
independent estimate of the influence of the small scales on the large scales it is necessary to 
estimate this quantity by some other means (possibly N-body simulations). 

7. Conclusions. 

We find that the lowest order nonlinear (~ e^) correction to the two point correlation 
calculated using the BBGKY hierarchy matches with the result obtained using the single stream 
approximation. Thus at this order there are no effects of multistreaming. The calculations 
presented here do not tell us if we we can get any effect due to multistreaming by going to higher 
orders of perturbation or whether the inability to follow the transition from a single streamed 
flow to a multistreamed flow is a limitation of the perturbative approach. In a separate study 
(Bharadwaj 1995) we have considered s system whose dynamics is governed be the Zel'dovich 
approximation instead of the gravitational dynamics. The transition from a single streamed flow 
to a multi-streamed flow occurs in the Zel'dovich approximation also. For this system we consider 
situations such that the initial flow is single streamed and it can be described by the gradient 
of a potential which is a Gaussian random field. We use two methods to perturbatively follow 
the evolution of the two point correlation function 1. The single streamed fluid equations which 
break down once multi-streaming occurs, and 2. distribution functions in phase space which are 
in principle valid even in the multi-streamed epoch. We find that the results obtained using 
these two methods match at all orders of perturbation indicating that there are no effects of 
multi-streaming at any order in a perturbative expansion. Non-perturbative analysis shows that 

A 

the effects of multi-streaming is of the form and hence they do not contribute in a 

perturbative expansion. We expect such conclusions to hold for the gravitational dynamics also. 
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It is known that if the initial power spectrum at small k has the form A;" with —3<n<—l, 
then the lowest order nonlinear correction to the power spectrum is a finite quantity which is the 
difference of two infinite quantities. This is inconvenient if one wants to numerically evaluate the 
correction to the power spectrum. We find that in real space, i.e. for the two point correlation 
function, this difficulty can be avoided by changeing the boundary conditions for the Laplacian 
operator. Usually when solving an equation of the form V^/(a;) = g{x) one uses the boundary 
condition limitx->oof(x) = 0. We find that if we use the boundary condition /(O) = then the 
above situation does not arise. 

To study the influence of small scales on the large scales we have considered two cases where 
the two point correlation is initially zero at large separations(i.e. $^^\x) = for x > xq). One of 
the cases has been treated analytically and the other numerically. We find that at large separations 
the induced two point correlation function has the form $^^\x) oc x~^. We also find that in both 
the cases the constant of proportionality is the same functional of the initial two point correlation 
function. 

We have numerically investigated some cases where the initial two point correlation has a 
power law form x~'^ at large x. The cases we have studied have a power spectrum of the form 
A;" with < < 3 at small k and have an exponential or gaussian cut off at large k. The cut off 
introduces a length scale and we find that at scales much smaller than this scale the nonlinear 
term enhances the growth of clustering. At intermediate scales we find that the nonlinear term 
changes sign more than once and it can act both ways i.e. to increase or decrease the clustering. 

At large scales we find that the behaviour of the nonlinear term depends on the condition 
whether $^^\x) oc $^^\x) or not. For the cases where this condition is satisfied we find that 
equation (78) gives a good fit to the lowest order nonlinear correction $^^\x) to the two point 
correlation. We find similar equations for the average of the two point correlation ^(^)(a;) and the 
pair velocity v^'^\x) too. We have interpreted equation (78) in terms of a simple diffusion process. 
For all the quantities the nonlinear term has the same sign as the corresponding linear term. 

For the case where $^^\x) is not proportional to $^^\x) we find that we obtain an equation 
similar to equation (78) with a different numerical coefficient. The quanitities $^'^\x) and v^'^\x) 
also have a similar behaviour, but for both of them it is not the leading linear part that contributes. 

In all the cases equation (78) shows the effect of the various scales on the large scales. 
Equation (78) has an integral of the initial two point correlation over all scales, and for the initial 
conditions that we have considered the small scales contribute the most. Thus we see that at the 
lowest order of nonlinearity the small scales effect the large scales. 

Hansel et al (1985) have studied the effect of large scales on the small scales and they find 
that in the weakly nonlinear regime the small scale perturbations get modulated by the large scale 
perturbations. This effect is like a diffusion in Fourier space because the effect of the disturbances 
at small k is to spread out the power spectrum at large k. It is very interesting that the effect of 
small scales on the large scales is a diffusion process in real space and the effect of the large scales 
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on the small scales is like a diffusion process in Fourier space. 

It is believed that perturbation theory is valid at a certain scale until ^(a;,f) at that scale 
is of the order of one. Another criterion for the applicability of perturbation theory is that the 
second order term should be smaller than the linear term. We have compared these two criteria 
and we find that the two break down at nearly the same epoch at small scales. At large scales we 
find that the second order term becomes of the same magnitude as the linear term at an epoch 
when ^(a;,f) <C 1. This happens because the second order term at large scales is influenced by the 
small scales. This happens at an epoch which is much later than the epoch when the perturbative 
treatment breaks down at small scales and because of the coupling of scales the second order term 
may be an overestimate of the actual nonlinear effects at large scales. For the CDM spectrum with 
the standard normalization this will occur only if we continue the evolution into the future. We 
also propose that for general initial conditions the true mean square displacement of the particles 
from their original positions rather than their linear value may be used to obtain a more realistic 
estimate of the influence of small scales on the large scales. 

We have tested the hypothesis that the dimensionless pair velocity h(^x,t) is an universal 
function of the average of the two point correlation ^(^x,t). In the linear regime we have 
h[x,t) = ^^[x,t). We find that the lowest order nonlinear correction introduces length scale 
dependent deviations from this relation and the pair velocity cannot be expressed as a universal 
function of the average of the two point correlation function. Although our investigation is valid 
only in the weakly nonlinear regime where the departures from the scaling relation are small, it 
indicates that as the evolution proceeds the contribution from the lowest order nonlinear terms 
will tend to enhance this deviation from the scaling relation. The approximate scaling relations 
found in N-body simulations by Hamilton et. al. (1991) seems to indicate that as the evolution 
proceeds, effects that have not been considered in our analysis tend to cancel out the deviations 
predicted by the lowest oreder nonlinear terms. This could be due to either higher order nonlinear 
effects or non-perturbative effects that come into play after multi-streaming. 

The author would like to thank Rajaram Nityananda for his advice, encouragement and for 
many illuminating discussions. 
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Fig. 1. — The linear two point correlation and its lowest order nonlinear correction for the case 
where the initial two point correlation correspnds to the self convolution of a sphere. The two 
curves have been given arbitrary displacements along the y axis for convenience of displaying. 
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Fig. 2. — The linear pair velocity and its lowest order nonlinear correction for the case where the 
initial two point correlation correspnds to the self convolution of a sphere. The two curves have 
been given arbitrary displacements along the y axis for convenience of displaying. 
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Fig. 3. — The initial two pint correlation ^^^\x) shown by the solid line and it average ^^^^{x) 
shown by the dotted line for two of the different initial power spectra considered.. The two plots 
on top show the small x behaviour. 
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Fig. 4. — The initial pair velocity for two of the different initial power spectra considered. The two 
plots on top show the small x behaviour. 
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Fig. 5. — The nonlinear correction to the two pint correlation $^^\x) shown by the solid line and 
it average ^(^)(a;) shown by the dotted line for two of the different initial power spectra considered. 
The two plots on top show the small x behaviour. 





Fig. 6. — The nonlinear correction to the pair velocity for two of the different initial power spectra 
considered. The two plots on top show the small x behaviour. 
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Fig. 8. — The solid line shows the two point correlation function ^{x,t) = ^^^\x,t) + ^^^\x,t) at 
cc = as a function of the scale factor S{t). The dotted line shows just the linear contribution. 
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Fig. 9. — The solid curves show the two point correlation function ^{x,t) = ^^^\x,t) + ^^^\x,t), 
its average ${x,t) and the pair velocity v(^x,t) at cc = .1 as a function of the scale factor S{t). The 
dotted line shows just the linear contirbution to these quantities. 
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Fig. 10. — The solidcurves show the two point correlation function ${x,t) = ^^^\x,t) + ^^^\x,t), 
its average ${x,t) and the pair velocity v(^x,t) at cc = 1.0 as a function of the scale factor S{t). The 
dotted line shows just the linear contirbution to these quantities. 





Fig. 11. — The solidcurves show the two point correlation function ${x,t) = ^^^\x,t) + ^^'^\x , t) ,its 
average ${x,t) and the pair velocity v(^x,t) at cc = 20. as a function of the scale factor S{t). The 
dotted line shows just the linear contirbution to these quantities. 
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Fig. 12.— The ratio for the initial power spectrum P{k) = e ^ is shown as a function of 

^{x,t) for different values of the scale factor S[t). S^^ corresponds to the value of the scale factor 
for which the nonlinear correction to the two point correlation at cc = is one tenth of the linear 
quantity {^'^^\x,t) ~ 0.1^'^^\x,t)). 
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Fig. 13.— The ratio for the initial power spectrum P(A;) = e is shown as a function of x for 

different values of the scale factor S{€). S^^ corresponds to the value of the scale factor for which 
the nonlinear correction to the two point correlation at cc = is one tenth of the linear quantity 



